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In this paper, we constrain the dimensionless Compton wavelength parameter Bq of f{R) gravity 
as well as the mass of sterile neutrino by using the cosmic microwave background observations, the 
baryon acoustic oscillation surveys, and the linear growth rate measurements. Since both the f{R) 
model and the sterile neutrino generally predict scale-dependent growth rates, we utilize the growth 
rate data measured in different wavenumber bins with the theoretical growth rate approximatively 
scale-independent in each bin. The employed growth rate data come from the peculiar velocity 
measurements at z = 0 in five wavenumber bins, and the redshift space distortions measurements 
at a = 0.25 and 2 = 0.37 in one wavenumber bin. By constraining the f{R) model alone, we get a 
tight 95% upper limit of logjg Bq < —4.1. This result is slightly weakened to logjg Bo < —3.8 (at 2(j 
level) once we simultaneously constrain the f{R) model and the sterile neutrino mass, due to the 
degeneracy between the parameters of the two. For the massive sterile neutrino parameters, we get 
the effective sterile neutrino mass < 0.62 eV (2cr) and the effective number of relativistic 

species < 3.90 (2ct) in the f{R) model. As a comparison, we also obtain < 0.56 eV 

(2a) and Neti < 3.92 (2a) in the standard ACDM model. 


One of the most profound puzzles in contemporary 
physics is the cause for the current cosmic acceleration 
[1, 2]. Although the mainstream for explaining this late¬ 
time cosmic acceleration is by introducing an extra dark 
energy fluid under the framework of general relativity (a 
typical dark energy candidate is cosmological constant A, 
corresponding to a ACDM cosmology) [3-5] , other possi¬ 
bility like the modification to general relativity on large 
scales also attracts more and more attention. A sim¬ 
ple modified gravity scenario is the so-called f{R) theory 
which generates acceleration via adding a function of the 
Ricci scalar R to the Einstein-Hilbert action [6, 7]. For 
recent reviews of f(R) theory, see, e.g.. Refs. [8-12]. 

Generally, f{R) theory introduces an extra scalar de¬ 
gree of freedom fn = df /dR characterized by a squared 
Compton wavelength proportional to fna = d^ f /dR?, or 
by a dimensionless quantity [13] 

Irr dR /' dlnH \~^ 

1 + /fl In a V ® ' 

where H = a/a is the Hubble parameter, a is the scale 
factor, and a dot denotes the derivative with respect to 
the cosmic time t. Due to the existence of this extra 
scalar degree of freedom, f(R) models tend to enhance 
the growth of the matter perturbations below the Comp¬ 
ton wavelength scale. Thus f(R) models generally pre¬ 
dict higher and scale-dependent growth rates compared 
with the standard ACDM model. By accurately measur¬ 
ing the growth rate of large-scale structure, we therefore 
have a powerful tool to constrain f(R) models. 

A direct way of measuring the large-scale growth rate 
comes from redshift-space distortions (RSD) in galaxy 
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surveys [14, 15], which depict the distortions of the mea¬ 
sured galaxy maps in redshift space. Such distortions 
arise from the peculiar velocities of galaxies, which ac¬ 
tually relate to the evolution of matter perturbation 
Sm- Thus RSD provide information about the large-scale 
growth rate /„ = dlnSm/dlna or the bias-independent 
growth rate fmO'siz) with as{z) the root-mean-square 
mass fluctuation in spheres with radius 8h~^ Mpc at red¬ 
shift z? Recent galaxy surveys have measured fmcrsiz) 
in different scale ranges with precision up to 7% [16-21]. 
Some works also have used these data to constrain some 
specific f(R) models [22, 23]. 

Nevertheless, current RSD models used to extract fm 
from the power spectrum of galaxy survey always assume 
a scale-independent fm on all scales, which is true for the 
standard ACDM model but not valid for the f(R) cos¬ 
mology. To see clearly, we plot the general fc-dependent 
growth rate fm{k, z) at 2 = 0 in Fig. 1 for the stan¬ 
dard ACDM and f(R) models with Bq = 10“"^, 10“® and 
10“®, where k is the comoving wavenumber and the sub¬ 
script “0” denotes the current value of the corresponding 
quantity. The plotting range for k covers most of cur¬ 
rent RSD surveys’ scales. It is clear that the growth rate 
in the f(R) model signihcantly deviates from the scale- 
independence. Particularly for the Bq = 10“'* case, the 
deviation arrives about 14% at A: = 0.3/iMpc“*, even 
larger than current RSD surveys’ precisions. So it will 


^ Note that, in practice, crgfz) is computed in linear regime. Al¬ 
though the expression, ag(z) = dkk'^P(k, z)W^(kr), 

where r = 8h~^ Mpc and P(k,z) is the matter power spectrum, 
contains an integral on all scales, the top-hat window function, 
W(kr) = 3(sin kr — kr cos kr) / (kr)^, decays rapidly when k is be¬ 
yond the linear regime, and thus the nonlinear modes contribute 
very little to the integral. 
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FIG. 1: The scale-dependent linear growth rate fm{k,z) at 
z = 0 for Bo = 10“^ (pink dashed), 10“® (red dot), and 
10~® (violet dashed dot) of f{R) model, and = 0-5 

eV (olive dashed dot dot) and 1.0 eV (blue short dashed) of 
sterile neutrino. The black solid line denotes the growth rate 
for the standard ACDM model. Note that we fix Qrn at the 
same values for all the curves, and A^efi = 4.046 for the curves 
of sterile neutrino. It is clear that ACDM model gives a scale- 
independent growth rate, while both f{R) and sterile neutrino 
predict scale-dependent growth rates. Besides, a larger Bq 
tends to enhance the growth rate more, while affects 

the growth rate in the opposite way. 


induce substantial systematic errors for the fit results 
when we simply use all the RSD data to constrain the 
f{R) models without careful selection according to their 
measured scales. Actually, by performing N-body sim¬ 
ulations, Ref. [24] illustrated that the measured growth 
rate from RSD surveys agrees with that predicted in f{R) 
model only for k < 0.06/iMpc“^ at z = 0. 

To avoid such systematic errors, we constrain /(i?) 
model using the scale-dependent growth rate data in 
this work. Recently, Ref. [25] tried to measure the 
scale-dependent growth rate by using the observations 
of peculiar velocities of galaxies from 6dF Galaxy Survey 
velocity sample in combination with a newly-compiled 
sample of low-redshift type la supernovae. The mea¬ 
surement obtained the scale-dependent growth rates in 
five k bins at z = 0: fmcrs{k) = 0.79 ± 0.21, 0.30 ± 
0.14, 0.32 ± 0.19, 0.64 ± 0.17, and 0.48 ± 0.22, for k e 
[0.005, 0.02], [0.02, 0.05], [0.05, 0.08], [0.08, 0.12], and 
[0.12, 0.15]/iMpc“^, respectively.^ Hereafter, we use PV 
to denote these five data points. Besides, we also use 
the growth rate measurements fm<^s{z) = 0.35 ± 0.06 


^ The original data are not gaussian. For the purpose of using 
these data to constrain f{R) model, we choose the upper limits 
of the data as the global errors. The growth rates are measured 
in a series of Ak ~ 0.03hMpc“^ bins in Ref. [25]. We assume 
that each data point is measured independently. We calculate 
the theoretical value of /mO’8(^i) each bin at ki = k^^^ with 
k^^^ the centre of the ith bin, and then compare it with the 
measurement value of the 2 th bin. 


at z = 0.25 and fmcrs{z) = 0.46 ± 0.04 at z = 0.37 from 
SDSS DR7 [19]. These two data points actually belong to 
the RSD measurement but are obtained from the power 
spectrum with k G [0.005, 0.033]hMpc“^. In this scale 
range, the growth rate in f(R) model can be approxima- 
tively considered scale-independent from Fig. 1. So it is 
appropriate to use these two data to constrain the f(R) 
model. Hereafter, we use RSD to denote these two data 
points. 

A successful f(R) model should give a stable acceler¬ 
ation solution at late times and also pass the local so¬ 
lar system tests. Several viable f(R) models that satisfy 
these conditions have been proposed in the literature [26- 
34]. However, we only care about the large-scale pertur¬ 
bations of f(R) models in this paper. So we consider a 
parametrized f(R) model whose background is the same 
as that of ACDM and perturbations are parametrized by 
two functions [35, 36], 


fx(k,a) 

l{k,a) 


1 1 -b |Af 

l-1.4-10-8|Ai|2a3 l-bA2 ’ 
1 + |A^ k'^a^ 

1 + |A^ k'^a^ ’ 


where the parameter Ai has dimension of length, which 
can be related to the current value of B via A^ = 
Bqc^/{2Hq) with c the speed of light. The functions 
li{k, a) and '-f{k, a) are used in the public MGCAMB code 
[37, 38] to quantify the modifications to the Poisson and 
anisotropy equations: 


k^'^ = —ir(fc, a)47rGa^{pA + 3{p + (3) 

/j2[<j) _ a)12TTGa^{p + P)cr, (4) 


where 41 and $ are two metric perturbation potentials, 
A is the comoving density perturbation, and cr denotes 
the anisotropic stress. It was showed that Eq. (2) to¬ 
gether with Eqs. (3) and (4) can well model the large 
scale growth in f{R) models for Bq < 1 [36]. In our pa¬ 
per, we will constrain Bq using above-mentioned growth 
rate data. Since the growth rate is very sensitive to the 
order of magnitude of Bq, we will take log,g Bq as a free 
parameter. 

Besides the value of Bq of f{R) model, another interest 
of this work is the mass of massive neutrinos. The total 
mass of massive active neutrinos in a f{R) universe has 
been widely reported in previous works, e.g.. Refs. [38- 
42]. So we here focus on the mass of massive sterile 
neutrinos. Recently, the extra massive sterile neutrino 
species has attracted much attention because of its magic 
power of simultaneously alleviating the tensions between 
different observations [43-52]. To describe the mass of 
sterile neutrino, we need two extra free parameters, the 
effective sterile neutrino mass and the effective 

number of relativistic species Nes. With these two pa¬ 
rameters, the true mass of a thermally-distributed sterile 
neutrino reads = (Aes - In 

order to avoid a negative N^g must be larger 
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than 3.046 in a universe with sterile neutrinos. The mass 
of sterile neutrino also corresponds to a characteristic 
scale, below which sterile neutrino is free-streaming and 
tends to suppress the growth of structure. Thus sterile 
neutrino can also lead to a scale-dependent growth rate 
(see Fig. 1). There might exist a degeneracy between Bq 
and the mass of sterile neutrino when we use the growth 
rate data to simultaneously constrain the two. We will 
constrain the parameters of sterile neutrino in the /(i?) 
model and verify whether this degeneracy exists or not. 

Our calculation is based on the pub¬ 
lic Markov-Chain Monte-Carlo package Cos- 
moMC [53]. The free parameter vector is: 

{Wfc, Wc, OmC, T, Us, ln(10l°As), logiQ^o, TO®®terile> 

where ujb, ujc, and Hq denote the background physical 
baryon density, physical cold dark matter density, and 
Hubble constant today, respectively, ^mc is approxima¬ 
tion to the angular size of the sound horizon at the time 
of last-scattering, t denotes the optical depth to reion¬ 
ization, and ln(10^*^Hs) and rig are the amplitude and 
the spectral index of the primordial scalar perturbation 
power spectrum. Following Ref. [40], we set a prior 
[—10, 1] for logiQ Bq. The priors of other free parameters 
are the same as those used by Planck Collaboration [54] . 
For massive active neutrino species, we assume its total 
mass = 0-06 eV. 

To constrain other cosmological parameters and break 
degeneracies between them, we also employ the observa¬ 
tional data from cosmic microwave background (CMB) 
and baryon acoustic oscillation (BAO) observations. For 
the CMB data, we use the temperature power spectrum 
data Cf^ [54] and lensing data [55] from Planck^ in 
combination with the WMAP 9-yr polarization (TE and 
EE) power spectrum data [56]. For the BAO data, we 
use the measurements from 6dFGS (z = 0.1) [57], SDSS 
DR7 (z = 0.35) [58], WiggleZ (z = 0.44, 0.60, and 0.73) 
[59], and BOSS DRll (z = 0.32 and 0.57) [60]. Note that 
the late-time integrated Sachs-Wolfe (ISW) effect and the 
lensing signal from the CMB measurements can also give 
constraint on Bq, since Bq affects the evolutions of the 
metric perturbation potentials $ and 4' from Eqs. (2)- 
(4), and the ISW effect and CMB lensing potential are 
directly related to 4> -|- 'k and $ -1- 4*, respectively. 

First, we constrain the f{R) parameter Bq without 
considering the extra sterile neutrino species. The main 
fit results are shown in Table I and Fig. 2. By using the 
CMB-I-BAO data and the scale-dependent growth rate 
PV data, we obtain the 95% upper limit of log^Q Bq < 
—2.7. Here note that the main contribution to the con¬ 
straint result of Bq comes from the PV data, since current 
CMB-I-BAO data (with the CMB lensing measurement) 


® After the completion of this work, the Planck 2015 results were 
reported (on arXiv). Note also that currently the Planck 2015 
likelihoods are not yet available. So the Planck data used in this 
work come from the Planck 2013 release. 
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FIG. 2: The one-dimensional posterior distributions 

and two-dimensional marginalized contours (68% and 95% 
CL) for the f{R) model from the CMB-I-BAO+PV and 
CMB-I-BAO+PV-I-RSD data combinations. 


can only give a weak constraint on Bq, e.g., Bq < 0.12 at 
2(7 level [61]. From Fig. 2, we can see a weak peak for 
log^o Bq at about —4 from the CMB-I-BAO-I-PV data. 
A peak with a nonzero Bq means that a relatively large 
growth rate is more favored by the data. The PV data in 
the first and fourth k bins may dominate the appearance 
of this peak, since the growth rates measured in these two 
bins are noticeably large. However, this peak for logj^g Bq 
disappears once we further add the RSD data into our 
analysis. Besides, the 95% upper limit of log^g Bq is im¬ 
proved to logiQ Bq < —4.1 with the help of the RSD data. 
These changes result from the fact that the RSD data fa¬ 
vor a lower growth rate and have higher precision than 
the PV data. 

Our fit result for Bq from the CMB-|-BAO-|-PV-(-RSD 
data is comparable to the limit, log^g Bq < —4.07 (at 
2(7 level), quoted in Ref. [40], where the WiggleZ galaxy 
power spectrum with /cmax = 0.2/iMpc“^ is used, and 
is also comparable to the limit from the recent Planck 
2015 result, log^g Hg < —4.01 (at 2a level), obtained 
by combining the RSD data at z = 0.57 [62]. These 
results represent the overall levels of constraint on the 
parametrized /(i?) model that can be achieved from the 
current information on linear scales. A tighter constraint 
on f{R) model can be obtained if the probes on smaller 
scales are utilized, such as the Planck Sunyaev-Zeldovich 
(PSZ) cluster counts [63], the strong lensing galaxy sig¬ 
nals [64], the distance indicators [65], and the fifth-force 
effects on diffuse dwarf galaxy components [66, 67]. How¬ 
ever, it is more difficult to control the possible systematic 
errors of these small-scale measurements. For example, 
the tension between the Planck temperature power spec¬ 
trum measurement and the PSZ cluster counts may be a 
hint of some unknown systematic errors existing in the 
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TABLE I: The mean values and la errors for the main parameters of the f(R) model, the ACDM+m®®tg,.iig + Ngg model, and 
the ALfi model from the CMB+BAO+PV and CMB+BAO+PV+RSD data combinations. For the parameters 

’^y*steriiei A^eff, and logjg Bq, we show their 95% upper limits and mean values in brackets. Note that the mass of sterile neutrino 
is in unit of eV and the Hubble constant Ho is in unit of km s“^ Mpc“^. 



f{R) 


+ + A^efF 

/(i?)+m®l,t5.rile + Neff 

Parameter CMB+BAO+PV 

+RSD 

CMB+BAO+PV 

+RSD 

CMB+BAO+PV 

+RSD 

sterile 

Neff 

logio Bo 

as 

Ho 

< -2.7(-6.1) 
0.3035 ± 0.0071 
0.850tg:°°« 

68.1 + 0.6 

< -4.1(-7.0) 
U.OUZ/ _Q 0075 

0.827in?8 
68.2 + 0.6 

< 0.43(0.18) 

< 3.96(3.47) 

0.3047 + 0.0080 

0.809lHg 

69.8lli 

< 0.56(0.23) 

< 3.92(3.43) 

U.JUU(_Q ggr^g 

0 700+0 031 
u. 1 »»_0.023 

69.611;? 

< 0.61(0.22) 

< 3.95(3.45) 

< -1.8(-5.6) 

u.auuu_Q QQgj^ 
n +0.032 

U.001_Q QY9 

69.8ll'i 

< 0.62(0.24) 

< 3.90(3.40) 

< -3.8(-6.8) 
0.3047l°o;“^i 

0.807l?;“l 

69.51];? 

In ^max 

4912.405 

4913.993 

4912.085 

4913.852 

4912.361 

4914.293 


PSZ measurement. Besides, the physics of f{R) theory 
on nonlinear scales is also complicated. So we do not in¬ 
corporate the information on small scales in our analysis. 

Next, we consider the cosmological constraints on 
f{R) parameter in a universe with sterile neutri¬ 
nos. The fit results from the CMB-fBAO-|-PV and 
CMB-t-BAO-t-PV-|-RSD data combinations are shown in 
Table I and Fig. 3. From Fig. 3, we also see a peak 
for logio ^0 at about —4 for the CMB-I-BAO+PV data. 
Comparing Fig. 3 with Fig. 2, we find that there indeed 
exists a slight degeneracy between the f{R) parameter 
and the sterile neutrino parameters. Due to this degen¬ 
eracy, the 95% upper limits of log^g Bq are enlarged to 
log^gRo < —1.8 and < —3.8 for the CMB-I-BAO-I-PV 
data and the CMB-|-BAO-fPV-l-RSD data, respectively. 
For the sterile neutrino parameters, we get < 

0.61 eV and N^g < 3.95 (2cr) from the CMB-kBAO+PV 
data, and w®®teriie < 0.62 eV and Ngg < 3.90 (2ct) 
from the CMB+BAO+PV-I-RSD data. To see how the 
f{R) gravity affects the fit value of the sterile neutrino 
mass, we also constrain the sterile neutrino parameters 
in the ACDM universe by using the same observational 
data. The main fit results are also shown in Table I. 
The CMB-I-BAO-I-PV data give w®®teriie < 0.43 eV and 
Aeff < 3.96 (2cr) in the ACDM universe. When the RSD 
data are also taken into account, the constraint results 
become < 0.56 eV and N^g < 3.92 {2a). We 

can find that the 95% upper limits of the effective ster¬ 
ile neutrino mass are also enlarged in the f{R) gravity, 
which also results from the degeneracy between the f{R) 
parameter and the sterile neutrino parameters. 

Current cosmic acceleration can result from either a 
dark energy fluid or a modification to gravity on large 
scales. Among various modified gravity theories, the 
so-called f{R) theory is a simple one. The key fea¬ 
ture of the f{R) theory that differs from the standard 
ACDM model is that f{R) models generally predict 
higher and scale-dependent growth rates. In this paper, 
we focus on a parametrized f{R) model, whose power 
of modification to the large-scale structure growth is de¬ 
scribed by a dimensionless Compton wavelength param¬ 
eter Bq. We constrain log^g Bq by using the CMB data. 
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FIG. 3: The one-dimensional posterior distributions and two- 
dimensional marginalized contours (68% and 95% CL) for the 
/(R)-|-m®®tg,.iig -I- Neff model from the CMB-I-BAO+PV and 
CMB+BAO+PV+RSD data combinations. 


the BAO data, and the linear growth rate data. To 
avoid the possible systematic errors when constraining 
the f{R) model from the growth rate data assumed to 
be scale-independent on all the survey scales, we utilize 
the scale-dependent growth rate data measured in sev¬ 
eral different wavenumber bins. They are the PV mea¬ 
surements of fm<y%{k,z = 0) in five wavenumber bins, 
k G [0.005,0.02], [0.02,0.05], [0.05,0.08], [0.08,0.12] and 
[0.12,0.150]/iMpc“^. We also use the RSD measure¬ 
ments of fmO'siz) at z = 0.25 and z = 0.37 in one 
wavenumber bin k G [0.005, 0.033]/iMpc“^. In each 
wavenumber bin, the theoretical growth rate predicted 
by the f{R) model can be approximatively considered 
scale-independent. By using the CMB+BAO+PV+RSD 
data, we get a tight 95% upper limit of log^g Bq < —4.1. 

We also extend our discussions to a universe with ster¬ 
ile neutrino species. Unlike the case in f{R) gravity, 
the extra sterile neutrino tends to suppress the growth 
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of structure below its free-streaming scale. We find 
a slight degeneracy between the f{R) parameter and 
the sterile neutrino parameters. Due to this degener¬ 
acy, the constraint on the f{R) parameter is slightly 
weakened to log^g Bq < —3.8 at 2a level from the 
CMB-I-BAO-I-PV-I-RSD data. For the massive sterile 
neutrino parameters, we get < 0.62 eV and 

Neff < 3.90 (2cr). As a comparison, we also obtain 
m®®teriie < 0.56 eV and A^eff < 3.92 (2a) in the stan¬ 
dard ACDM model. 
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